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1 Introduction 

One of the goals in a cosmic ray detection experiment is to determine the energy of the incident 
particles. Unfortunately, the primary energy cannot be measured directly at high energies where 
the flux is very low. Instead we take advantage of the cascades (or extensive air showers) produced 
by the cosmic rays in the atmosphere. The secondary particles that make up the cascade can be 
detected at ground level, as can the Cerenkov light that they produce. Alternatively, one can detect 
the atmospheric nitrogen fluorescence light induced by the passage of the shower. This technique, 
employed by the Fly's Eye detector and its successor the High Resolution Fly's Eye (HiRes) has 
the advantage that one can measure the number of charged particles in the shower as a function 
of depth in the atmosphere, N^hiX), where X is measured in g/cm?. We treat the atmosphere as 
a calorimeter, and the primary cosmic ray energy can be estimated by integrating the longitudinal 
profile Nch{X) and making corrections for "unseen" energy. 

In the past, the energy of a pure electromagnetic shower has been determined by Q, 

E f°° 

Eem = ^ Ne{X)dX, (1) 

where Xq is the electron radiation length in air, Ec is the critical energy of an electron in air, and A^e 
is the number of electrons in the shower. This equation implies that the electromagnetic energy is the 
total track length of all charged particles multiplied by an energy loss rate dE/dX given by Ec/Xq. 
One source of error here is the numerical value of the critical energy Ec, which has two definitions 
attributed to Rossi j^, and Berger and Selter Q. We turn to simulations to check the result and 
avoid this confusion. 

We do this using the modern shower simulation package CORSIKA Q. We simulate ultra-high 
energy showers complete with realistic fluctuations and realistic distributions of the energies of shower 



particles. An example is shown in Figure 1 where we plot the energy spectra of shower particles at 
the depth of the maximum size. Particle energies cover a wide range and lose energy to ionization at 
different rates (see inset to figure). 

For these simulations we replace Eq.(l) with a more general expression for the calorimetric energy, 



Ecai = a / N,h{X)dX, (2) 

JO 

where we integrate the charged particle longitudinal profile. We replace the constant in Eq.(l) with 
a parameter a representing the mean ionization loss rate over the entire shower. (This factor will be 
approximately equal to E^/ Xq and is calculated below). Given E^ai we must then make a correction 
to determine the cosmic ray energy Eq. The correction takes account of energy carried by high energy 
muons and neutrinos that ultimately deposit most of their energy in the ground. It also takes account 
of the small amount of energy that is lost to nuclear excitation. This "missing" energy has previously 
been parametrized by Linsley [0] and the Fly's Eye group |8|. 

In this paper we first describe some characteristics of the CORSIKA shower simulation package. 
We then use CORSIKA to simulate gamma-ray induced air showers (which have a very small "missing 
energy" component) to check the calorimetric energy method. Finally we simulate proton and iron 
induced showers to calculate the "missing energy" corrections for primary energies up to lO^'^eV. 

2 The Simulations 

CORSIKA is a versatile package for simulating air showers over a wide range of primary energies. 
Choices are available for the hadronic interaction model at the higest energies, and we have chosen 



the QGSJET |11] description which is in good agreement with Fly's Eye measurements. Within 



CORSIKA electromagnetic sub-showers are simulated with the EGS4 code. In EGS4, the cross 



sections and branching ratios are extended to lO^'^eV with the assumption that QED is vahd to these 
energies. In order to reduce CPU time, a thinning algorithm was selected within CORSIKA. That 
is, if the total energy of secondary particles from a given interaction falls below 10~^ of the primary 
energy, only one of the secondaries is followed, selected at random according to its energy Ei with a 
probability of pj = Ei/^jEj. The sum does not include neutrinos or particles with energies below 
the preset thresholds. In our simulations the threshold energies are 300, 700, 0.1 and 0.1 MeV for 
hadrons, muons, electrons and photons respectively. Particles below the threshold energies are not 
followed by the simulation. We chose an observation level 300 m above sea level and we simulated 
showers with zenith angles of 45 degrees. 

3 Calorimetric energy of an air shower 

Consider a purely electromagnetic shower. The primary particle energy E^^m can be approximated 

by 



AE{k)Me{k)dk, (3) 

where Afe{k) is the differential energy spectrum of electrons with kinetic energy k and AE(k) is the 
energy loss by each of those electrons in the calorimeter via ionization. This is only an approxi- 
mation because we have only included particles with kinetic energies above a threshold e. This is 
consistent with our simulations where we must impose a threshold of 0.1 MeV for photons, electrons 
and positrons. This integral can be carried out by summing over all the electrons produced in the 
simulation. 

We rearrange Eq.(3) and include the energy spectrum of particles as a function of atmospheric 
depth. 



K{k) = j^ Ar,(x)ne(fc,X)^^^ (4) 

where /\X{k) is the mean free path of electrons as a function of fc, N(,{X) is the total number of 
electrons at depth X, and ne(fc,X) is the normalized electron energy spectrum. Then, the electro- 
magnetic energy is approximated by 

Eem = j^ N,{X)i^J .^-{k)ne{k,X)dkjdX (5) 

The age parameter rather than the depth is often used to describe the stage of development of a 
shower. The energy spectrum of electrons can then be parametrized in terms of age. Since the age 
parameter is really only valid for a pure electromagnetic cascade, and since we will use the parameter 
in reference to hadronic showers, we will refer to our parameter as the pseudo age. We define it as 

gl^x) = 3 • (X - Xi) 

{X — Xi) + 2 • {Xmax — Xi) 

where Xi is the depth of first interaction and X^ax is the depth at which the shower reaches maximum 
size. Under this definition S{Xi) = 0, S{Xmax) = 1 and S{oo) = 3. 

One can then calculate the mean ionization loss rate {dE/dX) for the electrons in the shower 
(with energies > e) at age S", 

«(5) = J^^^{k)neik,S)dk (7) 

where Tie is now a function of 5. For comparison with Eq.(l) we rewrite Eq.(5) as: 



J-'em — 



<a>s Ne{X)dX, (8) 

JO 



where 



and < Nf, >ASi is the average number of electrons within a pseudo age bin ASj. 

We have simulated lO^^eV showers initiated by photons, protons and iron nuclei in order to 
calculate the mean energy loss rate over the entire shower, < a >s- We use bins of A^j = 0.1. 
Figure 2 shows a{S) as a function af age and we find < a >s is 2.186, 2.193 and 2.189 MeV/(g/cm^) 
for gamma, proton and iron induced showers respectively. All the errors in those < a >ss are less 
than 0.1%. 

This compares with the value of the ratio Ed Xq = 2.18 MeV/(g/cm^) used by the Fly's Eye 
analysis, where the values were taken to be Ec = 81MeV and Xo=37.1 g/cm^ [12]. This agreement 
may be a coincidence, since more recent values of the parameters from [|| are Ec = 86MeV (using 
Rossi's definition) and Xo=36.7 g/cm^, giving a ratio which is 7% higher than the typical simulation 
value of < a >5'. However, we note that the simulation results only include the energy loss rates for 
particles above the 0.1 MeV threshold. 

Figure 3 shows average shower profiles as a function of age for different primary masses and 
energies, with the shower size normalized to 1 at S* = 1. The difference in the average proton-induced 
shower profile at three different primary energies is smaller than the difference between the proton 
and iron average profiles at one energy. In other words, the shape of the shower development curve 
as a function of S is quite independent of primary energy or primary mass. It is also well known that 
for photon primaries the energy spectrum of the shower particles is a function of S only. We found 
it is also true for hadronic showers in our Monte Carlo study. Hence we can assume that our result 
for < a >5 can be applied over a range of primary masses and energies. 

We now apply Eq.(2) to some gamma-ray initiated CORSIKA showers, with a = 2.19MeV/(g/cm^). 
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The CORSIKA energy thresholds are set as described above with a threshold energy of 0.1 MeV for 
photons, electrons and positrons. We integrate the shower development curves in two ways for com- 
parison. In the first we numerically integrate the CORSIKA output which bins the development 
curve in 5 g/cvo? increments. Alternatively we fit a Gaisser-Hillas function (with variable Xq, X^ax, 
Xmax and A) to the CORSIKA output then integrate the function. Both methods give results which 
agree at a level of better than 1%. 

Table 1 shows the results for 500 showers. The calorimetric energy is about 10% lower than the true 
value. This is true even when we switch off processes which are not purely electromagnetic, namely 
/i+/i^ pair production and photo-nuclear reactions, which have small but important cross-sections 
in gamma-ray initiated showers. Two hundred such showers were generated and the results are also 
shown in Table 1, where we see that these showers have no muon content as expected. However the 
deficit in the calorimetric energy remains close to 10%. 

The solution to the problem is related to the simulation energy threshold of 0.1 MeV. We have made 
a detailed study with CORSIKA of the energy loss mechanisms and the characteristics of particles 
around 0.1 MeV. In particular we have summed the energy of particles that drop below the 0.1 MeV 
threshold. Table 2 shows that at 10"'^'^ eV, 88.4% of the primary energy is lost to the atmosphere 
through ionization by particles above 0.1 MeV. Electrons in the shower with energies below 0.1 MeV 
carry 9.0% of the primary energy, while sub-O.lMeV photons carry 1.2% of the primary energy. 
And the calorimetric energy derived by Eq.(2) is 88.8% of the primary energy, a good match to the 
ionization energy loss by particles above 0.1 MeV. 

We assume that the sub-0.1 MeV particles will eventually lose energy to ionization. The nitrogen 
fiuorescence efficiency is proportional to the ionization loss rate, so experiments like HiRes will detect 
light in proportion to the energy loss, even for very low energy particles. Thus the problem we 
experience with reconstructing the energy of CORSIKA simulations will not occur with the real 
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shower data. So for the further CORSIKA studies described below we have added 10% of the 
primary energy to the integrated energy loss result (from Eq.(2)) to take account of the sub-0.1 MeV 
particles that do not appear in the CORSIKA output. 

4 Energy estimation for hadronic showers 

We have described the calorimetric energy estimation for gamma-ray induced showers. We next 
consider hadronic showers where we expect the calorimetric energy to fall short of the primary 
energy because of so called "missing energy" - that energy channeled into neutrinos, high energy 
muons, and nuclear excitation. Much of this energy is deposited into the ground and is not visible 
in the atmospheric calorimeter. The first estimate of missing energy was obtained by Linsley who 
made measurements of electron and muon sizes at ground level and assessed the energy content of 
these components. The Fly's Eyes group parametrized Linsley's estimates as [Q]: 

Ecai/Eo = 0.990 - 0.0782 • E^^^-^"^^ (10) 

where Eq is the primary energy and E^-ai is the calorimetric energy derived from Eq.(2), both in units 
of 10^'^eV. This parametrization was said to be valid for 10^^eV< Eq < 10^'^eV. 

We have simulated proton and iron initiated showers at 8 primary energies from 3 x lO^^eV to 
10^°eV using CORSIKA. We apply Eq.(2) (with a mean energy loss rate of 2.19 MeV/ (g/cm^)) by 
fitting a Gaisser-Hillas profile to the CORSIKA development curve and extrapolating the profile to 
infinity and then integrating the function. We then add 10% of the primary energy to this result to 
take account of the CORSIKA threshold effect. Finally we compare this calorimetric energy with 
the primary energy, as shown in Figure 4. It is physically reasonable that the missing energy should 
decrease with increasing primary energy. Because of relativistic effects, charged pions produced in 
more energetic showers have an increased chance of interacting rather than suffering decay, reducing 



the fraction of energy immediately directed into muon and neutrino production. For comparison we 
also show Linsley's result. The solid line in Figure 4 shows the average behavior for proton and iron 
showers, which we express here as a function of Ecai (for practical convenience) in units of lO^'^eV, 

E^i/Eo = (0.959 ± 0.003) - (0.082 ± 0.003) • E~^f-^^°^°-°°^^ (11) 

which is valid for 3 x 10^^eV< Eq < lO^'^eV. Unfortunately it is never possible to know the primary 
particle mass on a shower-by-shower basis, so this average correction must be used. This lack of 
knowledge translates into an energy uncertainty, which is at most about 5% if the primary is hadronic. 
If the primary particle is a gamma-ray this assumption will overestimate the energy by up to 20%. Of 
course, if the shower development profile is obviously anomalous (as expected for gamma-ray showers 
above lO^^eV due to the LPM effect) gamma-ray primaries can be recognized and this systematic 
can be avoided. 

In an experiment like HiRes, atmospheric nitrogen fluorescence provides a measurement of ioniza- 
tion energy deposition, since the yield of fluorescence photons is proportional to this energy deposition 
Q. In the reconstruction process we convert the amount of light emitted by the shower at a particular 
depth to a number of charged particles, assuming that those charged particles are ionizing at the 
mean ionization rate which is a function of temperature and density |^] . This is taken into account in 
our analysis of real showers when we calculate the number of ionizing particles at a particular atmo- 
spheric depth. We then perform the path length integral (Eq.(2)), multiply by the mean ionization 
loss rate of 2.19MeV/(g/cm^)) and then make a correction for missing energy (Eq.ll). 

5 Conclusion 

We have re-investigated the veracity of estimating cosmic ray energy by using the atmosphere as 
a calorimeter. We have determined that, provided we use an appropriate mean energy loss rate. 



the technique provides a good estimate of primary energy for gamma-ray induced showers. For 
hadronic showers we have derived a correction function which accounts for energy not deposited 
in the atmosphere, so that the technique also returns a good estimate of primary energy for these 
showers. 
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With fi+fi- & -jN 


Without fi+fi' k jN 


So, eV 


Ecal/Eo 


^M 


N 


Ecal/Eo 


N^ 


N 


10" 


0.888 ± 0.004 


(2.146 ± 0.795)10^ 


(8.324 ± 0.392)10'' 


0.897 ± 0.003 


0.000 ± 0.000 


(8.448 ± 0.413)10*^ 


10" 


0.888 ± 0.005 


(2.823 ± 1.369)10* 


(7.881 ± 0.339)10^ 


0.898 ± 0.004 


0.000 ± 0.000 


(7.967 ± 0.392)10^ 


10^« 


0.889 ± 0.004 


(3.185 ± 0.916)10'^ 


(7.439 ± 0.271)10* 


0.898 ± 0.003 


0.000 ± 0.000 


(7.558 ± 0.281)10** 



Table 1: Results of CORSIKA simulations of gamma-ray induced air showers at three primary 
energies. The right-hand half of the table shows results from simulations where photo-nuclear and 
muon pair production processes have been switched off. The uncertainties shown are root mean 
squared errors. 



^0, eV 


Eioss/Eo 


^e(<0.1 MeV)/^o 


^^(<0.1 MeV)/^o 


Ecal/Eo 


1016 


0.888 ± 0.003 


0.090 ± 0.001 


0.010 ± 0.001 


0.888 ± 0.004 


10^^ 


0.884 ± 0.005 


0.090 ± 0.001 


0.012 ± 0.003 


0.888 ± 0.005 


IQis 


0.876 ± 0.007 


0.092 ± 0.002 


0.018 ± 0.005 


0.889 ± 0.004 



Table 2: Results from a study of energy conservation within CORSIKA. Gamma-ray showers were 
simulated at three primary energies Eq. Eioss refers to the energy lost to the atmosphere through 
ionization by charged particles with energies above 0.1 MeV. The fraction of the primary energy 
carried by sub-O.lMeF electrons and photons is shown in the next two columns. The fraction of 
primary energy determined by the calorimetric equation (final column) is consistent with Eiqss/Eq. 
Again, all uncertainties are r.m.s. 
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Figure 1: The mean energy spectra of photons, electrons and muons at 5 = 1 for 200 proton showers 
at 10^^ eV. The spike in the photon spectrum corresponds to electron-positron annihilation. The 
inset shows the energy loss rate (in MeV/g/cm^) by ionization of electrons in dry air over the same 
energy range as the main figure. The ESTAR code produced by US National Institute of Science and 
Technology (NIST) was used below 10 GeV|jT^ and this curve is extrapolated into the region above 
10 GeV. 
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Age,S 

Figure 2: The mean ionization loss rate dE/dX as function of S for gamma-ray, proton, and iron 
induced showers at lO^^eV. 
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Figure 3: The solid line is an average shower profile for 200 iron induced showers at 10^^ eV. The 
other three lines are average shower profiles for 200 proton showers at lO^^eV (short dashed line), 
lO^^eV (dotted line) and lO^^eV (long dashed line). 
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Figure 4: The functions for correcting the calorimetric energy to the primary energy, as a function 
of calorimetric energy. Shown are the corrections for proton showers (dotted hne) and iron showers 
(short dashed Une) and an average of the two(sohd Une). For comparison, Linsley's function is also 
shown. 
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